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Abstract. The TsaUis distribution has been used recently to fit the transverse 
momentum distributions of identified particles by the STAR [T] and PHENIX [5] 
cohaborations at the Relativistic Heavy Ion Collider and by the ALICE [3] and CMS [1] 
collaborations at the Large Hadron Collider. Theoretical issues are clarified concerning 
the thermodynamic consistency of the Tsallis distribution in the particular case of 
relativistic high energy quantum distributions. An improved form is proposed for 
describing the transverse momentum distribution and fits are presented together with 
estimates of the parameter q and the temperature T. 
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1. Introduction 

The Tsallis distribution has gained prominence recently in high energy physics with very 
high quality fits of the transverse momentum distributions made by the STAR [1] and 
PHENIX collaborations at the Relativistic Heavy Ion Collider and by the ALICE [3] 
and CMS |1] collaborations at the Large Hadron Collider. 

In the literature there exists more than one version of the Tsallis distribution [5l |6] 
and we investigate in this paper a version which we consider suited for describing results 
in high energy particle physics. Our main guiding criterium will be thermodynamic 
consistency which has not always been implemented correctly (see e.g. [3 El [9]). The 
explicit form which we use is: 
dN _ 
dpT dy 

(1) 



^ ^Pt^t cosh y 



, _ , niT cosh y — u 
l + (g-l)- 



T 

where pt and rriT are the transverse momentum and mass respectively, y is the rapidity, 
T and fi are the temperature and the chemical potential, V is the volume, g is the 
degeneracy factor. In the limit where the parameter q goes to 1 this reproduces the 
standard Boltzmann distribution: 
dN 

lim — = 

'j-S'i dpT dy 
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PtTUt cosh y ( uit cosh y - n\ 
[ T ) ■ 

In order to distinguish Eq. [1] from the form used by the STAR, PHENIX, ALICE and 
CMS collaborations [H El El H]. The motivation for preferring this form is presented 
in detail in the rest of this paper. The parametrization given in Eq. ([T|) is close (but 
different) from the one used by STAR, PHENIX, ALICE and CMS: 

d^iV diV (n-l)(n-2) / - niQ' 



dpTdy dy nC{nC + mo{n-2))\ nC ) ^' 

where n, C and mo are fit parameters. The analytic expression used in Refs. [H El |3l H] 
corresponds to identifying 

(4) 

g — 1 

and 

nC (5) 
q-1 

But differences do not allow for the above identifiction to be made complete due to an 
additional factor of the transverse mass on the right-hand side and a shift in the tansverse 
mass. They are close but not the same. In particular, no clear pattern emerges for the 
values of n and C while an interesting regularity is obtained for q and T as seen in 
Figs. and ([7]) shown at the end of this paper. 

In the next section we review the derivation of the Tsallis distribution by emphasizing 
the quantum statistical form and the thermodynamic consistency. 

2. Tsallis Distribution for Particle Multiplicities. 

In the following we discuss the Tsallis form of the Fermi-Dirac distribution proposed 
in [IDl [m [la Ha E] which uses 

^T^(^) = ^1^- (6) 

1 + exp^ [^) 

where the function expg(x) is defined as 

[1 + (g- if x>0 
[l + il-q)xf^^-''^ if x<0 

and, in the limit where g — )■ 1 reduces to the standard exponential: 



liniexpg(x) — )■ exp(x) 

The form given in Eq. (|6|) will be referred to as the Tsallis-FD distribution. The Bose- 
Einstein version (given below) will be referred to as the Tsallis-BE distribution [14] . 

All forms of the Tsallis distribution introduce a new parameter q. In practice 
this parameter is always close to 1, e.g. in the results obtained by the ALICE and 
CMS collaborations typical values for the parameter q can be obtained from fits to the 
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transverse momentum distribution for identified charged particles [3] and are in the 
range 1.1 to 1.2 (see below). The value of q should thus be considered as never being 
far from 1, deviating from it by 20% at most. An analysis of the composition of final 
state particles leads to a similar result [13] for the parameter q. 

The classical limit will be referred to as Tsallis-B distribution (the B stands for 
the fact that it reduces to the Boltzmann distribution in the limit where g — 1) and is 
given by result E] 



l + {q-'' 



' g-l 

T J ■ 

Note that we do not use the normalized g-probabilities which have been proposed 
in Ref. [6] since we use here mean occupation numbers which do not need to be 
normalized. In the limit where g — ?■ 1 all distributions coincide with the standard 
statistical distributions: 



limnf(E) =n'^(E), (9) 

g-i-l 

lim n^^ (E) =n^^{E), (10) 



\im n^^ (E) =n^^{E). (11) 

A derivation of the Tsallis distribution, based on the Boltzmann equation, has been 
given in Ref. [121117]. Numerically the difference between Eq. anfd the Fermi-Dirac 
distribution is small, as shown in Fig. ([1]) for a value of g = 1.1. 

The Tsallis-B distribution is always larger than the Boltzmann one if g > 1. Taking 
into account the large pt results for particle production we will only consider this 
possibility in this paper. As a consequence, in order to keep the particle yields the same, 
the Tsallis distribution always leads to smaller values of the freeze-out temperature for 
the same set of particle yields [15]. 

3. Thermodynamic Consistency 

The first and second laws of thermodynamics lead to the following two differential 
relations [18] 



de = Tds + fidn, (12) 
dP = sdT + ndfi. (13) 

where e = E/V, s = S/V and n = N/V are the energy, entropy and particel 
densities respectively. Thermodynamic consistency requires that the following relations 
be satisfied 

„ de 



ds 
de 



(14) 
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Figure 1. Comparison between tlie Fermi-Dirac and Tsallis-FD distributions as a 
function of the energy keeping the Tsallis parameter q fixed, for various values of 
the temperature T. The chemical potential is kept equal to one in all curves, the units 
are arbitrary. 



n 



dP 
djjL 
dP 
df 



(16) 
(17) 



The pressure, energy density and entropy density are all given by corresponding integrals 
over Tsallis distributions and the derivatives have to reproduce the corresponding 
physical quantities. For completeness, in the next section, we derive Tsallis 
thermodynamics using the maximal entropy principle and discuss quantum g-statistics 
in particular Bose-Einstein and Fermi-Dirac distribution by maximizing the entropy of 
the system for quantum distributions. This extends the derivation of Ref. [9]. We will 
show that the consistency conditions given above are indeed obeyed by the Tsallis-FD 
distribution. 
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4. Quantum Statistics 



= -gV 



The entropy in standard statistical meclianics for fermions is given in the large volume 
limit by: 

\ FD 1 FD 

- — -- n mn 

(27r)3 L 

+ (1 - n^^) ln(l - n^^)] , (18) 

where g is the degeneracy factor and V the volume of the system. For simplicity Eq. fll8l) 
refers to one particle species but can be easily generalized to many. In the limit where 
momenta are quantized, which is given by: 

S^"" = -9Y. K In ni + il- n,) ln(l - m)] , (19) 

i 

For convenience we will work with the discrete form in the rest of this section. The 
large volume limit can be recovered with the standard replacement: 

The generalization, using the Tsallis prescription, leads to [TOl [HI [12] 

St"" = -9ll K In, n, + {1 - n^f ln,(l - n,)] , (21) 

i 

where use has been made of the function 

x^~i — 1 

In^(^) ^ , (22) 

often referred to as q-logarithm. It can be easily shown that in the limit where the 
Tsallis parameter q tends to 1 one has: 

lim ln„(x) = ln(x). (23) 

g-s>l 

In a similar vein, the generalized form of the entropy for bosons is given by 

St"" = -9ll K K n,-{l + n,)« ln,(l + n,)] , (24) 

i 

In the limit g — )■ 1 equations [21] and [2H reduce to the standard Fermi-Dirac and Bose- 
Einstein distributions. Further, as we shall presently explain, the formulation of a 
variational principle in terms of the above equations allows to prove the validity of the 
general relations of thermodynamics. One of the relevant constraints is given by the 
average number of particles, 

j:n! = N. (25) 

i 

Likewise, the energy of the system gives a constraint, 

Y.<E, = E. (26) 
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it is necessary to have the power q on the left-hand side as no thermodynamic consistency 
would be achieved without it. The maximization of the entropic measure under the 
constraints Eqs. (125|) and (l26i) leads to the variational equation: 

5 



0, 



(27) 



where a and P are Lagrange multipliers associated, respectively, with the total number 
of particles and the total energy. Differentiating each expression in iEq. ([2] 



_5_ 

6ni 



6 



q-1 



rii 



-in , 



9-1 



n. 



q-l 



and 



9-1 



By substituting Eqs 



, ([2nD and ([3U]) into ([27D, we obtain 



9-1 



1 



[q-l 

Which can be rewritten as 



1 - nA'^'^ 



1 



-1 



1 — n,' 



Hi 



g-l 



l3Ei - a 



0. 



q-l 

and, by rearranging Eq. f l32|) . we get 

1 - Ui 



/3Ei + a, 



(28) 
(29) 

(30) 

(31) 
(32) 



[l + {q-l){/3E, + a)] — 



which gives the Tsallis-FD form referred to earlier in this paper as [TUl [HI [12] 

1 



Hi 



[l + {q-l){/3E, + a)]— + 1 
1 



exp (a + l3Ei 



+ 1 



(33) 



Using a similar approach one can also determine the tsallis-BE distribution by starting 
from the extremum of the entropy subject to the same two conditionss: 

6 



0, 



(34) 



which leads to 



1 



rii 



[l + (g-l)(/3E, + «)]—-!' 
1 



exp ((Ei - fi)/T) 



- 1 



(35) 



where the usual identifications a = -fJ^/T and P = 1/T have been made. 
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5. Proof of Thermodynamical Consistency 

In order to use the above expressions it lias to be sliown tliat tliey satisfy the 
thermodynamic consistency conditions. To show this in detail we use the first law 
of thermodynamics 

-E + TS + fiN 



P 



V 



(36) 



and take the partial derivative with respect to fi in order to check for thermodynamic 
consistency, it leads to 



dP 



1 

V 
1 
V 

+ 



_dE_ ^dS_ ^ dP[ 



T 

iv + V — — 



Ei — fi\ dni 
' T J 



then, by explicit calculation 



9/i T 
9^ " Y 



1 + (g - 1 



(37) 



T 



1-9 



1 + (g - 1 



T 



and 



'1 - Ui 



n. 



9-1 



1 + 



(g-l)(E,-/i) 



T 



Introducing this into equation [371 yields 
dP 



n, 



(3^ 



which proves the thermodynamical consistency [161 

We also calculate explicitly the relation in equation [TH can be rewritten as 



dE 
'dS 



n dT 



dE I dE dfi 
_ dT dfi dT 
~ dS_ _i_ as dp ' 

dT ' dfi dT 

since n is kept fixed one has the additional constraint 
dn dn 

= ar^^ + djl"^^ = °' 

leading to 



(39) 
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Now, we rewrite [39] and HQ] in terms of the following expressions 



and 



dE 

&r 

dE 
dS_ 

&r 

dS_ 

dn 



dT V 



dn 1 



Y.QE,n. 



q- 



-i duj 
dT'' 

-i dnj 




dfi 



By introducing the above relations into equation 
becomes 

dE dEdf^_^ ^_^dn^ 
&T^ df^dT ^^'"'''^ dT 



(fEj {ni 



nj) 



dnj dn^ 



9/i dT 



Y^qn] 



Y.^E,{n,n, 



,9-1 



-^iduj 
dfi 



.iduj 



Where the abbreviation 



Cij = [niTij 



>9-l 



9/i 

dni duj duj dui 



the numerator of equation 



(41) 



dT dfx dT dfx 

has been introduced. One can rewrite the denominator part of equation [39] as 



(42) 



dS dS dfx 
&T ^ 'dJldT 



n'j dj 



^iduj 



djj 



q J2iEi - /i) (ninjY ^ dj 
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where 



-nf-' + (1 - n,Y-^ (E, - /i) 



-n 



q-l 



q-1 rp . 

hence, by substituting equation HH and US] in to [321 we find 

EiCij 

T- 



dE 

as 



since ^ Cij = 0, this finally leads to the desired result 



(44) 



dE 

as 



(45) 



Hence thermodynamic consistency is satisfied. 

It has thus been shown that the definitions of temperature and pressure within the 
Tsallis formalism for non-extensive thermostatistics lead to expressions which satisfy 
consistency with the first and second laws of thermodynamics. 

6. Tsallis Fit Details 

The total number of particles is given by the integral version of Eq. fl25|l . 



N = gV 



(2^ 



1 + (g - 1 



E- 



19/(1-9) 



T 



(46) 



The corresponding (invariant) momentum distribution deduced from the equation above 
is given by 



r.dN 1 

E = qVE 

(Pp ^ (27r)3 



(47) 



which, in terms of the rapidity and transverse mass variables, becomes 



dN 



dyprdpT 



gv 



tut cosh y 
(27r)2 

1 + (g - 1) 



ruT cosh y — II 
T 



9/(1-9) 



(48) 



At mid-rapidity y = and for zero chemical potential this reduces to the following 
expression 



dN 



dpT dy 



gV 



y=0 



Pttht 



1 + (g - 1 



19/(1-9) 



(49) 



Fits using the above expressions based on the Tsallis-B distribution to experimental 
measurements published by the CMS collaboration [3] are shown in Figs. (2), (3) and (4) 
are comparable with those shown by the CMS collaboration. Fits to the results obtained 
by the ALICE collaboration [3] are shown in Fig. (5) for iT'',K~,p. The resulting 
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parameters are considerably different from those obtained from Eq. [3] and are collected 
in Table I. The most striking feature is that the values of the parameter q are fairly 
stable in the range 1.1 to 1.2 for all particles considered at 0.9 TeV. The temperature 
T cannot be determined very accurately for all hadrons but they are consistent with a 
value around 70 MeV. 

For clarity we show these results also in Fig. (|6]) for the values of the parameter q and 



Particle 


1 


T (GeV) 


7r+ 


1.154 ±0.036 


0.0682 ±0.0026 




1.146 ±0.036 


0.0704 ± 0.0027 


K+ 


1.158 ±0.142 


0.0690 ±0.0223 


K~ 


1.157 ±0.139 


0.0681 ± 0.0217 




1.134 ±0.079 


0.0923 ±0.0139 


P 


1.107 ±0.147 


0.0730 ± 0.0425 


P 


1.106 ±0.158 


0.0764 ±0.0464 


A 


1.114 ±0.047 


0.0698 ± 0.0148 




1.110 ±0.218 


0.0440 ± 0.0752 



Table 1. Fitted values of the T and q parameters for strange particles measured by 
the ALICE [3] and CMS collaborations ^ using the Tsallis-B form for the momentum 
distribution. 

in Fig. d?]) for the values of the Tsallis parameter T. The striking feature is that the 
values of q are consistently between 1.1 and 1.2 for all species of hadrons. 

7. Conclusions 

In this paper we have presented a detailed derivation of the quantum form of the Tsallis 
distribution and considered in detail the thermodynamic consistency of the resulting 
distribution. It was emphasized that an additional power of q is needed to achieve 
consistency with the laws of thermodynamics [H]. The resulting distribution, called 
Tsallis-B, was compared with recent measurements from the ALICE [3J and CMS 
collaborations |1] and good agreement was obtained. The resulting parameter q which 
is a measure for the deviation from a standard Boltzmann distribution was found to 
be around 1.11. The resulting values of the temperature are also consistent within the 
errors and lead to a value of around 70 MeV. 
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(GeV/c) 

Figure 3. Comparison between tlie measured transverse momentum distribution for 
A as measured by the CMS collaboration 4 and the Tsallis-B distribution. The full 
line is a fit using the parameterization given in Eq. (j49p to the 0.9 TeV data with the 
parameters listed in Table I. 
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Figure 4. Comparison between the measured transverse momentum distribution for 
E!~ as measured by the CMS cohaboration [4, and the Tsalhs-B distribution. The fuU 
Una is a fit using the parameterization given in Eq. (j49p to the 0.9 TeV data with the 
parameters listed in Table I. 
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Figure 5. Comparison between the measured transverse momentum distribution 
for vr^, and p as measured by the ALICE collaboration [3] and the Tsallis-B 
distribution. The lines are fits using the parameterization given in Eq. to the 0.9 
TeV data with the parameters listed in Table I. Full line is for tt^, the dotted line is 
for , the dashed line is for anti-protons. 



The Tsallis Distribution in Proton-Proton Collisions at \/s = 0.9 TeV at the LHC. 15 
1.5 



1.4 - 
1.3 - 
1.2 
1.1 
cr 1 
0.9 
0.8 
0.7 
0.6 
0.5 



5 5 



p-p 900 GeV 



J. 



J. 



J. 



7t 



J. 



A E" 



K K s 
Figure 6. Values of the Tsallis parameter q for different species of hadrons. 
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Figure 7. Values of the Tsallis temperature T for different species of hadrons. 



